Local Potential Model of the Hoyle Band in 12 C 
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We describe the excited + state of 12 C at 7.654 MeV, often called the Hoyle state, in terms of a 
local potential 8 Be+a cluster model. We use a previously published prescription for the cluster- 
core potential to solve the Schrodinger equation to obtain wave functions for this state, and also for 
higher angular momentum states of the same system. We calculate energies, widths and charge radii 
for the resulting band of states, with particular emphasis on the recently discovered 2 + state. We 
examine various choices of the oscillator quantum number 2n + L = G for the cluster-core relative 
motion, and find that G = 6 leads to the most coherent description of the properties of the states 
and is consistent with recent experimental data on the L = 2 state. 

PACS numbers: PACS index numbers: 21.10.-k, 21.60.6x, 21.10.Tg and 27.20.+n 



I. INTRODUCTION 

When considering how carbon could be produced in 
stellar nucleosynthesis, Hoyle |l] famously predicted that 
the 12 C nucleus must have an excited + state in the 
vicinity of the 3-a breakup threshold. Such a state was 
duly found [2] in short order on the basis of Hoyle's sug- 
gestion. Its existence was essential for 12 C to be pro- 
duced at an adequate rate to open up the pathways to 
the synthesis of still heavier nuclei |3j . Not only had this 
state escaped both experimental detection and theoreti- 
cal prediction up to that point, but it has continued to 
pose severe challenges to nuclear structure models ever 
since. 

The shell model struggles to describe any low-lying ex- 
cited + states in p-shell and sd-shell nuclei. To this day, 
the most advanced no-core shell model has not succeeded 
in reproducing the excitation energy of the Hoyle state 
[U • Models invoking a 3a chain 0-0] were able to accom- 
modate the excitation energy, but could not reproduce 
the decay width of the Hoyle state. Better success was 
enjoyed by orthogonality condition model calculations of 
the 3a system, using a semi-microscopic a — a interac- 
tion 0, IH , which gave a suitably sized a-decay width in 
the 8 Be(0+)a channel. This suggested that the domi- 
nant structure of the Hoyle state must be 8 Be(0^)a in 
a relative s-state. Such a conclusion was backed up by 
fully microscopic 3a calculations using the Resonating 
Group and Generator Coordinate Methods |10l4l3l | . More 
recently fully microscopic calculations using Antisym- 
metrized Molecular Dynamics (AMD) [l4| and Fermionic 
Molecular dynamics (FMD) [15| have been able to give 
a good account of the low- lying spectrum of 12 C without 
assuming alpha clustering a priori. For the particular 
case of the very loosely bound Hoyle state a three-a con- 
densate wave function [l6| has also been shown to have a 
large overlap with the FMD wave function. The current 
situation is well summarized in a recent review pTj . 

In parallel with these theoretical developments, new 
life has recently been breathed into the experimental pro- 



gram. Given that the Hoyle state has a large mean square 
radius, it is natural to suggest the possibility of rotating 
this structure and producing a band of excited states with 
a large moment of inertia. Although Friedrich et al. [l8[ 
claim there are no such rotational states, most 3a models 
do support them, with the 2 + state generally expected 
to have an excitation energy in the region of 10 MeV. 
Initial searches for this 2 + state via beta decay [l9| were 
discouraging, but in recent times three separate experi- 
ments have provided evidence of its existence. Itoh et al. 
[2(j see a 2+ state at 9.9 MeV, with a width of 1.0 MeV, 
via inelastic a-particle scattering in their 12 C(a, a') mea- 
surements. Freer et al. see a 2 + state at 9.6 MeV, 
with a width of 0.6 MeV, in their inelastic proton scatter- 
ing 12 C(p,p') work. Gai [22] confirms the existence of a 
2 + state around 10 MeV (but without being able to mea- 
sure a width) in photonuclear disintegration 12 C(7, 3a) 
studies. 

In view of this resurgent interest in an excited Hoyle 
band we ask in this paper to what extent the known 
data can be accounted for (and further excited states pre- 
dicted) by a local potential 8 Be+a cluster model. This 
model provides a physically transparent and calculation- 
ally straightforward description of the energies, widths 
and charge radii of the known states in the proposed 
Hoyle band. It also throws some light on the question of 
how high in angular momentum such a band might con- 
tinue. In this way it can be a useful guide to experimental 
groups searching for as yet unidentified higher L states. 
It is also illuminating to see how far one can get with a 
simple but physically motivated model. Although ideally 
the full rigours of a more microscopic and computation- 
ally intensive approach might seem preferable, a rather 
simplified effective nucleon-nucleon potential is needed to 
carry this to fruition. It might be that a phenomenologi- 
cal approach can produce a band of states with properties 
closer to the experimental values. In any event it is in- 
teresting to compare the results from different theoretical 
treatments with each other as well as with experiment. 

In the next section we describe our local potential 
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8 Be+a cluster model. After that, we compare the 
model's results with the available data and discuss the 
possible existence of additional states. Finally, we sum- 
marize our conclusions. 



II. LOCAL POTENTIAL CLUSTER MODEL 

The cluster model employed here was first proposed 
to study the excited 4p-4h band of 16 O (bandhead at 
6.05 MeV) and the ground state band of 20 Ne [H] and 
has subsequently been applied to a wide range of nuclei 
across the Periodic Table from 6 Li to 242 Cm. However, it 
has not previously been used for 12 C, because the ground 
state band of that nucleus is known to be oblate (from 
the sign of its quadrupole moment [24]]) and a two-body 
cluster-core system inevitably produces a prolate defor- 
mation. However, the indications that the Hoyle state, 
and any associated rotations of it, have a large 8 Be(0^)a 
component, marks out the Hoyle band as ideal territory 
for our local potential two-body cluster model. 

In general, we model a nucleus as two even-even sub- 
nuclei of mass A x and A 2 separated by a distance R, in- 
teracting through a deep, local nuclear potential Vn(R) 
and a Coulomb potential Vc{R) appropriate to a point 
cluster and a uniformly charged spherical core. The nu- 
clear part has previously been parametrised in the form 



V N (R) = -V 



[l + exp((i?-i? )/a))] 
1-x ) 



[l + cxp((i?-i? )/3a)] 3 



(1) 



with parameter values given by 



v = 54.0 MeV, a = 0.73 fm, and x = 0.33 (2) 



and Vq related to vq by 



V = 



A X A 2 



{A! + A 2 ) 



Vq. 



(3) 



The radius parameter Ro of the potential is determined 
by fitting to the experimental energy of the spectrum's 
band head, and is linked to the choice of G = 2n + L 
(see below). This potential is now incorporated into the 
cluster-core relative motion Hamiltonian Hq (R) , and the 
resulting Schrodinger equation 



i? (R)$GnL(R) = E GnL $GnL(B.) 

UGnL(R) 



E, 



GnL " 



R 



Y LM (6,<j>). (4) 



is solved numerically. We label the wave functions and 
energies with the global quantum number G = 2n + L, 
where n is the number of internal nodes in the radial 
wave function and L the orbital angular momentum. 



We must choose the value of G large enough to guar- 
antee that the Pauli exclusion principle is satisfied by 
excluding the cluster constituents from states occupied 
by the core nucleons. For the 8 Be + a system this re- 
quires G > 4. The appropriate choice of G is not clear cut 
when we are not using an oscillator potential, although 
we can certainly use oscillator considerations as a guide. 
The lowest allowed value, G = 4, would correspond to 
packing the a nucleons into the p-shell, which seems an 
unlikely description of an excited state in 12 C. Neverthe- 
less, for completeness, we present calculations employing 
G = 4,6 and 8, with a view to deciding a posteriori on 
the basis of a comparison with experimental data which 
is the best choice. 

Solving the Schrodinger equation produces excitation 
energies and their associated widths directly. The result- 
ing wave functions can also be used to calculate mean 
square charge radii of the states from 



(Z 1 + Z 2 )(R 2 ( 12 C)) = Z 1 (R 2 ( s Be)) + Z 2 (R 2 (a)) 



(Z X A\ + Z 2 A\) 
(Ax + A 2 ) 2 



(5) 



where (Z X ,A X ) = (4,8) and (Z 2 ,A 2 ) = (2,4) for the 8 Be 
+ a system. To apply this formula for charge radii we 
need to supplement the wave functions generated above 
with a description of the ground state of 8 Be (or at least, 
of its mean square charge radius). An excellent descrip- 
tion of a — a scattering phase shifts and 8 Be has already 
been given within the local potential cluster model [261 ] 
using a Gaussian nuclear potential 



V N (R) = V G exp(-aR 2 



(6) 



with V G = 122.6225 MeV and a = 0.22 fm" 2 and a 
Coulomb potential 



V C (R) 



z x z 2 cTi(m) 

R 



(7) 



with (3 — 0.75 fm~ and where Z x = Z 2 = 2 for the 
a — a system. We adopt this description of 8 Be wholesale 
because, within our model, it can hardly be improved 
upon. 

As a consistency check on the widths of the states we 
can make use of our earlier work (see for example [27j ) on 
charged particle decay widths. Within the two-body local 
potential model, a semiclassical approximation leads to 
an a-partial width of 



T Q = F — exp -2 / dR k(R) , (8) 
V v Jr * J 

where R 2 and i?3 are the two outermost classical turning 
points and fi is the reduced mass of the system. The 
normalization factor F is given to good accuracy by 



F 



! dR 
Rl 2^) 



1, 



(9) 
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with i?i the innermost turning point, and the wave num- 
ber k(R) is 

k(R)=(j£\Q-V(R)\) 12 (10) 

and Q is the experimental energy of the decaying state 
relative to the two-body breakup threshold. 

III. ENERGIES, WIDTHS AND RADII OF 
HOYLE BAND STATES 

As outlined in the previous section, we calculate 8 Be 
+ a cluster states, using a previously published prescrip- 
tion for the potential [25|], without adjusting any of the 
free parameters. We identify the + state with the Hoyle 
state, and go on to calculate similar states of higher an- 
gular momenta. Each band of states is labelled by the 
choice of G = 2n + L. The maximum possible value of 
L in this scheme is equal to G itself (corresponding to 
the nodeless wave function) . We determine the potential 
radius i?o of the potential by fitting so as to reproduce 
the experimental energy of the + Hoyle state exactly. 
The resulting values for i?o are listed in Table I. It is 
also useful to know the position and maximum height of 
the potential (Coulomb barrier) since this gives an upper 
limit on the energy of a resonant state with the corre- 
sponding value of L. Therefore we also list these values 
in Table I for the three values of G = 4, 6 and 8 under 
consideration. 



TABLE I: Potential Maxima for G = 4, 6 and 8 





Vmax 


(MeV) at i? max 


(fm) 


L 


G = 4 


G = 6 


G = 8 




Ro = 0.8883 fm R = 1.9386 fm R a = 2.9008 fm 





1.63 at 6.20 


1.42 at 7.25 


1.26 at 8.28 


2 


3.12 at 5.16 


2.44 at 6.41 


2.02 at 7.56 


4 


9.23 at 3.48 


5.74 at 5.28 


4.24 at 6.62 


6 




13.87 at 4.02 


8.82 at 5.74 


8 






17.26 at 4.86 



Figure 1 shows the potentials (nuclear, Coulomb and 
centrifugal combined) obtained from the procedure de- 
scribed above for G = 4, 6 and 8. It is clear from in- 
spection that it would be no surprise to discover that the 
state with the highest expected L- value for each value of 
G was completely unbound (or at best precariously res- 
onant). In fact, none of the states of the Hoyle band is 
bound. They are all (at best) resonances, and in some 
cases rather close to the top of the Coulomb barrier, so 
the calculation of their energies and widths is numerically 
delicate. In view of this we have found it expedient to 
cross check our results using a variety of different calcu- 
lational methods. 



• We use our own bound state code to perform a 
numerical integration of the Schrodinger equation, 
but rounding off the potential at its maximum value 
for radii in excess of that radius where the maxi- 
mum is achieved (i.e. V — V max for R > R max ). 
This serves to give a good estimate of the energies 
which can be used to provide a starting energy for 
the methods described below, and also as a consis- 
tency check on these subsequent results. 

• We use the published code GAMOW [Hj], 
which employs complex arithmetic to solve the 
Schrodinger equation and fits Gamow tails to the 
resonant states, so as to evaluate energies and 
widths of the states. We need this code princi- 
pally for the wave functions it generates which can 
be used eventually in the calculation of the mean 
square charge radii. 

• We use the published elastic scattering code SCAT2 
(29j which allows us to monitor the scattering phase 
shifts of the 8 Be + a system as a function of the 
centre of mass (CM) energy, and thereby to identify 
the peak energies and widths of the resonant states 
in the various partial waves. 

• As an order of magnitude check, we also evalu- 
ate the widths of the states using the semiclassical 
method discussed in the previous section, Eq.(8). 
Previous experience using this approach to calcu- 
late half-lives for charged particle decay suggests 
that it can be expected to agree with the true value 
to within about a factor of two. 

The code SCAT2 writes out the scattering matrix ele- 
ments rji in the form 1 — Re(?7;) and Im(ry;), where the rji 
are related to the transmission coefficients T\ by 

T, = l-\ m \ 2 . (11) 

We run the code at successively incremented CM ener- 
gies, centred on the values indicated by our own bound 
state code and by GAMOW, and monitor the behaviour 
ot the scattering matrix as the energy increases through 
the suspected resonance region. As the energy of a 
resonant state is approached from below, the value of 
1 — Re (77;) for the appropriate partial wave I rises from 
near zero through 1 towards 2. We take the energy range 
over which it rises from 0.5 to 1.5 as the width of the 
resonant state. The value of Im(?7i) also rises from 
towards 1, and then falls again to zero. By systemat- 
ically working through the energy regions indicated by 
the bound state code (and taking successively finer incre- 
mental energy grids as necessary), we obtain the width 
values reported for the L > states in Table II. The 
widths for the + states were obtained using the semi- 
classical aaproximation described in the previous section, 
Eqs.(8-10). We note that the experimental width of the 
+ state is 8.5 ± 1.0 eV, with which the corresponding 
theoretical value obtained with G — 6 in Table II is in 
good agreement. 
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The energies from all these methods are mutually com- 
patible and we present the average of them as the peak 
energy of each resonance in Table II. We have added 7.365 
MeV to the CM scattering energies to obtain excitation 
energies relative to the 12 C ground state. (Note that the 
+ state's energy of 7.654 MeV was used to fit the po- 
tential radius Rq for all three values of G) 



TABLE II: Calculated state energies, widths and (R 2 ) for 
G = 4, 6 and 8 



code and confirmed by GAMOW using the potential of 
Ref.[26[ in the formula 



L 


G = 4 


E(MeV) ± T 
G = 6 


G = 8 





7.654 ± 3(eV) 


7.654 ± 6(eV) 


7.654 ± 13(eV) 


2 


No state 


9.61 ± 360(keV) 


8.85 ± 68(keV) 


4 


No state 


13.71 ± 1.24(MeV) 11.52 ± 300(keV) 


6 




No state 


16.25 ± 660(keV) 


8 






No state 


(r 2 ) 


11.688 fm 2 


13.457 fm 2 


15.553 fm 2 




3.42 fm 


3.67 fm 


3.94 fm 



Figure 2 shows the energy dependence of 1 — Ke(rji) 
and Im(rji) for the G — 6,L = 2 resonance near the CM 
energy of about 2.22 MeV. This is both a typical case, 
and the state that we are most interested in (because 
of recent experimental results [20|-[22j which have at last 
found a rotational state of the Hoyle band) . The best de- 
scription of this state within our model is achieved using 
G = 6. For G = 4 there is no resonant state at all, and for 
G = 8 the excitation energy is too low and the width too 
narrow. We are not unduly concerned that this value of 
G does not tally "nicely" with oscillator shell model con- 
siderations, which might suggest G=8 since most of the 
low-lying intruder states in this region are described as 
4p-4h excitations in the shell model. We are not using a 
harmonic oscillator potential and will generate rather dif- 
ferent wave functions with exponential rather than Gaus- 
sian tails, and so do not expect an identity of G values 
between the two cases. Our calculations also predict an 
excited 4 + state at 13.71 MeV. We predict a width of 
1.24 MeV for it, which would make it hard to detect in 
experiments. We note however that this excitation en- 
ergy places the state above the barrier maximum for an 
L = 4 state in a G = 6 band, indicating that the result 
is not completely reliable, although still indicative of a 
possible 4 + state in this region. 

Previous analyses of inelastic electron scattering from 
12 C [3(} indicate that the Hoyle state has an abnormally 
large charge radius. We therefore calculate the values of 
(r 2 ) implied by our model for G = 4,6 and 8. We first 
need to calculate a mean square charge radius for 8 Be. 
We do this by using the mean square charge radius for an 
a-particle of 1.6757 fm 3l| and the mean square separa- 
tion of the two a-particles calculated by our bound state 



(R 2 )* Be = (R 2 ) a + -(R 2 



(12) 



This yields a mean square charge radius for 8 Be of 10.634 
fm 2 (with a square root of 3.261 fm). This, in turn, serves 
as input for the mean square charge radius of the Hoyle 
state in the formula 



2 12 

{R )i2C-Hoylc = )a + -^(R }bc 

+ —{R 2 )a-Bc 

y 



(13) 



Our G = 6 value of 3.78 fm is somewhat below the value 
of 3.87 fm deduced by Funaki et al. [3(3] but significantly 
larger than the values deduced from earlier GCM calcu- 
lations of 3.50 fm [l2[ and 3.47 fm [ll| obtained from full 
3a calculations using Volkov effective two-nucleon forces. 



IV. CONCLUSIONS 

It has long been suspected that the Hoyle state in 12 C 
might have rotational excitations built upon it so that a 
Hoyle band could be present in the 12 C spectrum. Re- 
cent experimental work has located a 2 + state a little 
below 10 MeV with a width of about 600 keV which is 
a strong candidate for such a structure. This has moti- 
vated us to apply a local potential 8 Be-a cluster model, 
with previously published potential prescription, to the 
system to see if we can reproduce the existing data and 
predict properties of other similar states. 

The calculation is numerically delicate but obtaining 
closely similar excitation energies from three different 
methods give us a good degree of confidence in our re- 
sults. We find that with a global quantum number of 
G = 6 we are able to give a good account of the width 
and root mean square charge radius of the Hoyle state it- 
self and a reasonable description of the excitation energy 
and width of the proposed excited 2 + member of the pu- 
tative Hoyle band. We note that our calculated 2 + state 
is very close to the top of the Coulomb barrier. We also 
predict a rather wide 4 + state at an excitation energy of 
roughly 14 MeV, and find that the band certainly termi- 
nates here (if not already with the 2 + state). If the 4 + 
state is truly 1.25 MeV wide, it will be almost impossi- 
ble to see against the background of a known 4 + state 
at 14.083 MeV with a width of 258 keV 32]. However, 
our calculation is not completely reliable for this state, 
and there may be a case to be made that the known 
4 + state at 14.083 MeV should be assigned to the Hoyle 
band. Perhaps there are even two 4 + states there. This 
is an interesting conundrum which further experimental 
investigation might be able to resolve. 
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FIG. 1: Local potentials (nuclear + Coulomb + centrifugal) for the 8 Be + a system corresponding to all possible L-values for 
the global quantum numbers G = 4, 6 and 8. The potentials for G — 4 are displaced upwards by 100 MeV and those for G = 8 
downwards by 100 MeV to aid visibility. 
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FIG. 2: The energy dependence of 1 — Re(rn) and Im(r?j) for the G = 6, L = 2 resonance near the CM energy of 2.22 MeV 
obtained from the code SCAT2. The energy scale has been shifted slightly to align it with the average energy calculated for 
the 2 + excitation by our different calculations. 



